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Abstract: We analyze the dynamical behavior of a quasi-isotropic Universe in 
CO ■ the presence of a cosmological fluid endowed with bulk viscosity. We express the 

viscosity coefficient as a power- law of the fluid energy density: C — Co e s ■ Then 

■ we fix s = ] /2 as the only case in which viscosity plays a significant role in the 
\ singularity physics but does not dominate the Universe dynamics (as requested 

— , ■ by its microscopic perturbative origin). The parameter Co is left free to define the 

intensity of the viscous effects. 

■ Following the spirit of the work by E.M. Lifshitz and I.M. Khalatnikov on the 
[~^. \ quasi-isotropic solution, we analyze both Einstein and hydrodynamic equations 

up to first and second order in time. As a result, we get a power-law solution 
I existing only in correspondence to a restricted domain of Co- 

' 

H ! 1 Introduction 

i 

Near the cosmological singularity, the isotropic nature of the Universe corresponds to a class 
of solutions of the Einstein equations containing three physically arbitrary functions of the 
space coordinates. In the case of a radiation-dominated Universe, such a class was found by 
E.M. Lifshitz and I.M. Khalatnikov in 1963 I and then generalized to an arbitrary fluid state 
equation in [2] . Earlier extensions of this quasi-isotropic scheme were provided in j3j [4J, [5l [6] , 
where different evolutionary stages of the Universe are characterized. 

The original work by Lifshitz and Khalatnikov treated the quasi-isotropic model as a Tay- 
lor expansion of the 3-metric tensor in powers of the synchronous time. However, further 
investigation outlined the necessity of treating generic power-law components of the 3-metric. 
We fix our attention on the relevance of dealing with bulk viscous properties of the cosmolog- 
ical fluid approaching the Big-Bang singularity. In fact, bulk viscosity is a phenomenological 
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issue inherent to the difficulty for a thermodynamical system to follow the equilibrium config- 
uration. Since asymptotically near the Big-Bang the volume of the Universe has a very fast 
time variation, we naturally expect the bulk viscous effects arise in the dynamics. It is worth 
noting that such viscous contributions can not dominate the Universe evolution because of 
their thermodynamical perturbative origin; nevertheless, we are interested in those regimes 
where such effects are not at all negligible. 

General analyses of the Universe behavior in presence of bulk viscosity were given in 
[51 H] in which the bulk viscosity coefficient C is described by a power-law of the energy 
density of the fluid, i.e., C = Co £ s (Co 7 s = const). As far as this phenomenological ansatz is 
referred to the early isotropic Universe, it is easy to realize that the choice s = 1 /2 prevents 
dominating viscous effects. On the other hand, simple considerations, as well as the analysis 
presented in the works by J.D. Barrow [TOl HTJ [TJ], indicate that the case s < 1/2 leads to 
negligible contributions of the viscosity to the asymptotic regime towards the Big-Bang. As a 
consequence, in studying the singularity physics, the most appropriate form of the power-law 

is C = Co Ve- 
in [13], the stability of a flat Friedmann - Robertson - Walker (FRW) Universe was in- 
vestigated including bulk viscosity in the cosmological fluid dynamics (see also [2]). As a 
result, it was shown that the behavior of the perturbations is damped, forward in time, by 
these new features with respect to the standard case of an ideal fluid. In the present work 
we generalize the investigation to the whole class of quasi-isotropic solutions. Our aim is to 
determine the conditions on the parameter £0 (i.e., on the viscosity intensity) which allows 
for the existence of a quasi-isotropic regime for the radiation dominated Universe. 

For this purpose, we investigate the Einstein equations under the assumptions proper of the 
quasi-isotropic model. We separate zeroth- and first-order terms into the 3-metric tensor and 
the whole analysis follows this scheme of approximation. In the search for a self-consistent 
solution, we make use of the hydrodynamics equations in view of fixing the form of the energy 
density. As a result, we prove the existence of a quasi-isotropic solution, which has a structure 
analogous that provided by Khalatnikov, Kamenshchik and Starobinsky [SJ. Of course, in our 
solution the power-law for the leading 3-metric term is sensitive to the viscosity parameter 
Co- In particular, we find that such a solution exists only if when £0 remains smaller than 
a certain critical value. Finally, we determine the density contrast and its dependence on 
Co- This behavior confirms and generalizes the result obtained in |13j about the damping of 
density perturbations by the viscous correction. 

The scheme of the paper is as follows. In Sec. 2, a review of the quasi-isotropic solution [T] 
is presented; after setting the form of the 3-metric expansion, the integration of the Einstein 
equations is performed in order to obtain a solution for the energy density, the density 
contrasts and the 3-velocity of the perfect fluid filling the space time. In Sec. 3, we write 
down the expression of the 3-metric expansion in a more general form without fixing the time 
power-law exponents. The construction of the generalized Ricci tensor follows. In Sec. 4, we 
introduce the viscous-matter contribution into the Einstein equations and then the energy- 
momentum tensor is written down in presence of bulk viscosity. Furthermore, the Taylor 
expansion of the energy density is addressed using the 00-components of the gravitational 
equations. The viscous parameter s is here set equal to 1 /2 and a brief analysis of its physical 
meaning is performed. In Sec. 5, some comments on the validity of the addressed model are 
presented. We firstly discuss the role of the microphysical horizon concerning the thermal 
equilibrium of the expanding Universe and the related hydrodynamical treatment. Then we 
analyze in some details the effects of the first- (shear) viscosity and show how it must be 
neglected in the dynamics, for the appearance of internal consistences in the model. Finally, 
a comment on the introduction of a causal thermodynamics into the very early dynamics is 
addressed. In particular, we treat the viscous fluids description by using a relaxation- equation. 
In Sec. 6, we integrate the 00-Einstein equation matched together with the hydrodynamical 
one. The integration is faced separately for the zeroth- and first-order analysis and a critical 
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threshold for the viscous parameter £o is found by imposing the constraint for the consistence 
of the adopted model. As a result, we get the expressions of the energy density and the density 
contrast in presence of bulk viscosity. In Sec. 7, we integrate the Oa- and a/3-gravitational 
equations and an expression of the 3-velocity as a function of the perturbed metric term is 
found. Brief concluding remarks follow. 

2 Review on the LK-quasi-isotropic solution 

In 1963, E.M. Lifshitz and I.M. Khalatnikov pQ first proposed the so-called quasi-isotropic 
solution. This model is based on the idea that the space contracts maintaining linear-distance 
changes with the same time dependence order by order (i.e., a Taylor expansion of the 3- 
metric is addressed). In this approach, the Friedmann solution becomes a particular case of 
a larger class of solutions existing only for space filled with matter [TS] . 

The metric evolution is strongly characterized by the matter equation of state. For an ultra- 
relativistic perfect fluid, characterized by the relation p — e/3, the spatial metric assumes the 
form "f a f3 ~ a a /3 t, asymptotically as t — > (the cosmological singularity is set by convention 
in t = 0), where a a g are assigned functions of the coordinates. As a function of time, 
the 3-metric is expandable in powers of t. The quasi-isotropic solution is formulated in a 
synchronous system (i.e., go a — 0, goo = —1), which is not strictly a co-moving one. The 
line element writes as 

ds 2 = -dt 2 + j af) (t, x" l )dx a dx (:S , (1) 
with a spatial metric of the form 

7a/3 = ta a0 + t 2 b aP + ... , 7 a/3 = t- 1 a alB - b af3 , (2) 

where a a/3 is defined as a a/3 ap.y — 6"; furthermore, the relation 6^ = a" 7 6 7/ 3 is ensured by 
the scheme of approximation. 

The Einstein equations in the synchronous system assume the form [16] 

= 3 + I K a K p = — \T , (3a) 

Ri = 5^ W7 *g) + PS = Tg ~\TSi, (3c) 
where c=G=l and the tensor n a p and its contractions read 

n af 3 = da al 3 = a a(3 + 2 1 b a p , n p a = 7 /3<5 n a s = t' 1 6% + , k = d t In ^7 = 3 r 1 + b , 

(4) 

and 7 = det (70-/3) ~ t 3 (l + tb) det(a a p). Matter is described by an ultra-relativistic perfect 
fluid energy-momentum tensor 

T^v = (p + e) u^u v +pg liV = | {Au^ + g^) , (5) 

which provides the following identities 

T ° = \ e (~Au 2 + 1) , T° = fe UQU °, l£ = - § e t^u" , T = 0. (6) 

Calculating the left-hand side of (j3"a|) , (|3b[) up to zeroth- 0(l/t 2 ) and first-order 0(1/ 1) in 
l/t, we rewrite them respectively as 

3 b e 2 1 ,8 s 4e 

~4^ + 2l = 3 ( ° + ); 2 ( ;Q ~ = — 3- ( 7 ) 
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(8) 



Because of the identity — 1 = u^u^ ~ — Uq + t~ x u a ug a a ^ , it is immediate to see that e ~ t~ 2 
and u a ~ i 2 ; hence, in the asymptotic limit t — > 0, itg ~ 1 (uo = — 1). From the first equation 
of (|Tp. one can find the first two terms of the energy density expansion, while, from the second 
equation, the leading term of the velocity arises 

- JL A - * 2 

Because of ([5]) , the expression for the density contrast S can be found as first and zeroth-order 
energy density ratio, i.e., 

5 = -jbt. (9) 

This behavior implies that, as expected in the cosmological standard model, the zeroth-order 
term of energy density diverges more rapidly than the perturbations and the singularity is 
naturally approached with a vanishing density contrast in this scenario. 

Besides the solutions for e and u a , one has to consider the pure spatial components of the 
gravitational equation (|3c|) . Up to first approximation, the Ricci tensor can be written as 
P@ = P£/t, where P@ is constructed by the constant 3-tensor a a g. The terms of order t~ 2 
automatically cancel out, while those proportional to t^ 1 give 

H + \K + &&*2 = o. (io) 

Performing the trace of this equation, a relation between the quite arbitrary six functions 
a a g and the coefficients b a g from the next-to-leading term of expansion can be determined: 
= — A /3pP + 5 /i8 P5P. It is worth reminding that, in the asymptotic limit t — > 0, the 
matter distribution becomes homogeneous because e approaches a value independent of b. 

Now, using the Ricci identity P^.g — 1 /2p : p, the useful relation 6^.^ = 7 /9b :a can be 
determined; this gives the final expression for the 3-velocity distribution as 

t 2 

u a = — b- a . (11) 
y 

This result implies that, in this approximation, the 3-velocity is a gradient field of a scalar 
function fixed by the non perturbed metric a a p. As a consequence, the curl of the velocity 
vanishes and no rotations take place into the fluid. 

Finally, it must be observed that metric (|2|| allows for an arbitrary 3-space coordinate 
transformation and the solution above contains only 6 — 3 = 3 arbitrary space functions 
arising from a a g. A particular choice of this functions, those which correspond to the space 
of constant curvature {PP ~ 6@), can reproduce the pure isotropic and homogeneous model. 



3 Generalized quasi-isotropic line element 

In order to generalize the quasi-isotropic solution of the Einstein equations for the presence 
of dissipative effects into the evolution of the energy source, we deal with a more complex 
(no longer in integer powers) form of the 3-metric ^ HI [J. In this respect, we take the 
spatial metric as 

l a( i = t x a aB + i" b a3 , ^ = t~ x a aP - i»- 2 * b a ? . (12) 

Here, the constraints for the space contraction (i.e., x > 0), and for the consistence of the 
perturbation scheme (i.e., y > x) have to be imposed for the proper development of the 
model. In this approach, the extrinsic curvature and its contractions read 

K a3 = xt x ~ x a a3 + yt y ~ 1 b a g , (13a) 

ni=xt~ 1 5P + (y-x)ty~*- 1 bi, (13b) 

K =3xi- 1 +(y-x)t y - x - 1 b, (13c) 
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furthermore, we calculate the following useful relation 



8t = \ k = | xr 1 + \ (y - x) t v - x ~ l b . (14) 

We are now able to write down the final form of the Ricci-tensor components contained in 
the Einstein equations (J3J). These new expressions allow us to generalize the original quasi- 
isotropic approach. Our aim is to obtain constraints and relations for the exponents x, y in 
order to guarantee the existence of the solutions of our model. They explicitly read 

*8 = + (y - *) & - > ( 15a ) 

_ x(3x-2) {y-x){2y + x-2) (y-x)x p fg Pf_ 

H oc- U 2 °* + tfl-y+x b <* + U 2-y+x b6 a + p + t 2x-y ' ^ bC > 

We note that in equation ()15cp . P@ represents the 3-dimensional Ricci tensor constructed 
by the metric a a p. On the other hand, the higher-order term P*@ denotes the part of Pg 
containing the 3-tensor b a p. 



4 The form of energy density in the viscous approach 

In the quasi-isotropic solution the Universe is assumed, in correspondence with the standard 
cosmological model, to be described by the energy-momentum tensor of an ultra-relativistic 
perfect fluid. In connection with the development of new cosmological models, the discovery 
of the cosmic acceleration suggests matter to play an essential role at different stages of 
cosmological evolution and it can obey very different equations of state |17) . Thus, corrections 
in this sense to the original formulation of the quasi-isotropic solution can be useful in this 
new context. 

In this work, we treat the immediate generalization of LK scheme. We consider the pres- 
ence of dissipative processes within the fluid dynamics, as expected at temperatures above 
O(10 16 GeU); this extension is represented by an additional term in the expression of the 
energy-momentum tensor ([5]) and it can be derived from thermodynamical properties of the 
fluid [HEiiilQl]. The new tensor reads 

T^u = (P + e)u,j,u v +pg» v = | (4^^ + g^) -Cuf p (Ufj,u v + g^ v ) , (16) 

P = p-(u? pl (17) 

where p = e/3 denotes the usual thermostatic pressure in correspondence of an ultra- 
relativistic state equation and £ is the bulk viscosity coefficient. In this work, we neglect 
the so-called shear viscosity (first viscosity) since, in the case of quasi-isotropic cosmological 
evolution, the displacement of the matter layers with respect to each other appears only in 
a higher-order analysis (this kind of viscosity represents the energy dissipation due to such 
an effect). 

The coefficient £ has to be expressed in terms of the thermodynamical parameters of the 
fluid. In particular, here we assume this quantity as a function of the Universe energy density; 
according to literature developments 110[ lllj , we express C as a power-law of the form 

C = Coe s , (18) 

where Co is a constant and s is a dimensionless parameter whose behavior was discussed by 
V. A. Belinskii et al. [3 E] for asymptotic values of the energy density. This study yields the 
constraint ^ s ^ 1 /i for its range of variation in correspondence of large values of e. 
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Let us now write the expressions of the mixed components of tensor (|T6|) up to higher-order 
corrections as 

T ° = -f ( 4u o - 1) + Coe s «f M («g - 1) , (19a) 

T =-3( e s u^ , (19b) 

Tf = | (4^ + «£) - Coe s + Si) , (19c) 

T° = | e u Q u° - Coe s u& u Q u° , (19d) 
where the divergence of the 4-velocity reads as 

= d t In ^7 = f a* -1 + § (V ~ x) ty-^ 1 b . (20) 

Here, we assume, as in the non- viscous case, the relation: Uq ~ 1 (with uo = —1), whose con- 
sistence must be verified a posteriori comparing the time behavior of the quantities involved 
in the model. Taking into account expressions ()19a|) (Il9b|) . we can recast now the Einstein 
equation (|3a|) in the form 

3x(2 — x) , . , „ , 6 9i > . 3(y — x) M .„ . 

~ 4,2 + ^-^)(^- 1 )^2-^ = - £ + ' T7h-7^' (21) 

In what follows, we fix the value s = 1 / , 2 in order to deal with the maximum effect that 
bulk viscosity can have without dominating the dynamics. In fact, the notion of this kind of 
viscosity corresponds to a phenomcnological issue of perturbations to the thermodynamical 
equilibrium [5], [T5] . In this sense, we remark that, if s > 1 /2, the dissipative effects become 
dominant and non-perturbative. Moreover, if we assume the viscous parameter s < 1 /2, the 
dynamics of the early Universe is characterized by an expansion via a power-law a(t) ~ i 2 / 37 
starting from a perfect fluid Friedmann singularity at t = (here 7 is identify by the relation 
p = (7 — 1) e). After this first stage of evolution, where viscosity does not affect at all the 
dynamics, the Universe inflates in the limit t — > 00 (i.e., out of our approximation scheme) to 
a viscous deSitter solution characterized by alt) ~ e Hot , H = ^o/3 = l /z (CoV / 3/7) 1 ^ 1_2s - 1 

snug. 

Since, in this work, we deal with the asymptotic limit t — * 0, we only treat the case 
s = 1 /2 in order to quantitatively include dissipative effects in the primordial dynamics. 
From equation (j2~l|) . if s = 1 /2, we expand the energy density e as follows 

e ei6 r ^ f ei b \ 

where the constants eo, ei are to be determined combining the 00-gravitational equation with 
the hydrodynamical ones, comparing the terms order by order, as treated below. We remark 
that only for the case s = 1 /2 all terms coming on the left-hand and the right-hand side 
respectively of equation (|2ip result to have the same time behavior up to first order because 
of ~ 



5 Comments on the adopted paradigm 

In this Section, we discuss in some details the hypotheses at the ground of our analysis of 
the quasi-isotropic viscous Universe dynamics. In particular, we investigated some peculiar 
features of the very early evolution (near the cosmological singularity) since their presence 
leads to a specific treatment of the viscous phenomena. 

It is well known |20j the crucial role played in cosmology by the microphysical horizon, 
as far as the thermodynamical equilibrium is concerned. In the isotropic Universe, such a 
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quantity is fixed by the inverse of the expansion rate, H^ 1 = (a/a) (a being the scale factor 
of the Universe and the dot identifies time derivatives) and it gives the characteristic scale 
below which the elementary-particle interactions are able to preserve the thermal equilibrium 
of the system. Therefore, if the mean free-path of particles £ is greater than the microphysical 
horizon (i.e., £ > H^ 1 ), no real notion of thermal equilibrium can be recovered at the micro- 
causal scale. If we indicate by n the number density of particles and by a the averaged 
cross section of interactions, then the mean free-path of the ultra-relativistic cosmological 
fluid (in the early Universe the particles velocity is very close to speed of light) takes the 
form £ ~ 1/na. Interactions mediated by massless gauge bosons are characterized by the 
cross section a ~ a 2 T~ 2 (g — y/4iva being the gauge coupling strength) and the physical 
estimation n ~ T 3 leads to the result £ ~ l/a 2 T [20] . During the radiation-dominated era 
H ~ T 2 /m P i, so that 

l~-£—H- x . (23) 

Therefore in the case of T > a 2 mpi ~ O(10 16 GeV), i.e., during the earliest epoch of pre- 
inflating Universe, the interactions above are effectively "frozen out" and they are not able to 
maintain or to establish thermal equilibrium. To complete this consideration we remark that, 
at temperatures grater than O(10 16 GeU), the contributions to the estimation above due to 
the mass term of the gauge bosons can be ruled out for all known and proposed perturbative 
interactions. 

As a consequence of this non-equilibrium configuration of the causal regions characterizing 
the early Universe, most of the well-established results about the kinetic theory [211 EH E3] 
concerning the cosmological fluid nearby equilibrium become not applicable. Indeed all these 
analysis are based on the assumption to deal with a finite mean free-path of the particles 
and, in particular, results about the characterization of viscosity are established when pure 
collisions among particles are retained. However, when the mean free-path is grater than the 
micro-causal horizon (which, in the pre-inflating Universe, coincides with the cosmological 
horizon), £ can be taken of infinite magnitude for any physical purpose. 

The fundamental analysis of the viscous cosmology is due to the Landau school pj [HI [H] ; 
since they were aware of these difficulties for a consistent kinetic theory, such an analysis 
was essentially based on an hydrodynamical approach. A real notion of the hydrodynamical 
description can be provided by assuming that an arbitrary state is adequately specified by 
the particle- flow vector and the energy momentum tensor alone [24]. In particular, the en- 
tropy flux has to be expressible as a function of these two hydrodynamical variables without 
additional parameters. Following this point of view, the viscosity effects are treated on the 
ground of a thermodynamical description of the fluid, i.e., the viscosity coefficients are fixed 
by the macroscopic parameters which govern the system evolution. In this respect, the most 
natural choice is to take such a (shear and bulk) viscosity coefficients as power-laws in the 
energy density of the fluid (for a detailed discussions see [7])- Such a phenomenological as- 
sumption can be reconciled, for some simple cases, with a relativistic kinetic theory approach 
|25j . especially in the limits of small and large energy densities. 

Addressing the hydrodynamical approach, we are lead to retain the same equation of 
state which would characterize the corresponding ideal fluid. This fact is supported by idea 
that the viscosity effects provide only small corrections to the thermodynamical setting of the 
system. As clarified above, in the present analysis, we deal with the case in which bulk viscous 
corrections are of the same order of the perfect fluid contributions, in order to maximize their 
influence in the Universe dynamics. Nevertheless, since we are treating an ultra-relativistic 
thermodynamical system, which is very weakly interacting on the micro-causal scale, there 
are well-grounded reasons to describe it by the equation of state p = e/3. 

Another important point concerning the ground assumptions of our model, is why the shear 
viscosity (rj) is not addressed in the present scheme. Indeed, this kind of viscosity accounts 



7 



for the friction forces acting between different portions of viscous fluid. Therefore, as far as 
the isotropic character of the Universe is retained, the shear viscosity must not provide any 
contributions, as discussed in [5]. On the contrary, the rapid expansion of the early Universe 
suggests that an important contribution comes out from the bulk viscosity as an averaged 
effect of a quasi-equilibrium evolution. 

Indeed, our present analysis deals with small inhomogeneous corrections to the background 
FRW-metric. Thus, at first-order in our solution, shear viscosity should be, in principle, 
included into the dynamics. In this sense it is shown in [7] that, if the bulk viscosity coefficient 
behaves like £ ~ e s , then the corresponding shear one behaves as r\ ~ e p , where p must satisfy 
the constraint condition p ^ s + 1 /2. Here we treat the case s — thus getting p ^ 1 for the 
rj coefficient. This issue is incompatible with the symmetries and the approximations here 
addressed. In fact, the shear viscosity provides, among others, an equivalent contribution to 
the bulk one, since the energy-momentum tensor of the viscous fluid contains the term 

T/jo, ~ ... - (C - %rj) u? p (u,j,u v + g pv ) + ... . (24) 

We now observe that, at zeroth-order, u? p ~ 0(l/t), while the first-order correction in 
the energy density behaves like 0(l/t x ) and we will show the relation 1 ^ x < 2 in Sec. 
6. Since the request x ^ 1 comes out from the zeroth-order analysis, which by isotropy 
is independent of the shear contribution, we can conclude that, for our model, the shear 
viscosity would produce the inconsistency associated to the term 0(l/t px+1 ). The point is 
that the request px + 1 2 would make such a contribution dominant in the model, against 
the basic assumption. Thus, to include shear viscosity in a quasi-isotropic model, we should 
choose the case s < 1 /2 which is out of the aim of this paper since it is devoted to maximize 
the bulk effects in a coherent cosmological dynamics. 

To conclude this section, we would like to discuss the question concerning the implemen- 
tation of a causal thermodynamics for our cosmological model. Indeed, the hydrodynamical 
theory of a viscous fluid is applicable only when the spatial and temporal derivatives of the 
velocity of the matter are small [26j [27l [24] . This condition is necessarily violated in the 
asymptotic limit near the cosmological singularity. This way viscous fluids would have to be 
described by using a relaxation equation similar to the Maxwell equation in the theory of 
viscoelasticity [5]. 

In this scheme, the energy-momentum tensor assumes the form 

Tpv = e u p u u - (p + a) (g^ + u p u u ) , (25) 

where p denotes the thermostatic pressure and a is the bulk-stress density. In the very early 
Universe, the relation between a and the relaxation time To reads as follow 

a + ar - C< P ■ (26) 

The relaxation time can be expressed as tq/C ~ l/ e • this physical assumption follows from 
the fact that the transverse-wave velocity in matter has finite (nonzero) magnitude in the 
case of large values of e [9] . 

In this scheme, we are able to express the time dependence of To. Since, at leading order, 
e ~ 1/t 2 , we obtain, using Eq. (fT5|) . the following behavior for the relaxation time To ~ i 2 ~ 2s . 
In our model we deal with the case s = 1 /% which yields To ~ t and, if we address a power-law 
dependence on a (according the structure of the solution) such as & ~ a/t [§], relation (T2H|) 
rewrites as 

cT = (oe s u? p . (27) 

From this analysis we can apply the standard expression for the bulk viscous hydrodynamic 
taking into account the reparameterization £o — > Co of the bulk coefficient. 

The considerations above allow us to regard the subtle paradigm of the causal thermody- 
namics, having in mind that it would affect only qualitative details of our analysis, but it 
could not alter the validity of our results. 
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6 Solutions of the 00-Einstein and the hydrodynamical 
equations 

In Section 4, we exploited Eq. (|3a|) in order to obtain the qualitative expression for the energy 
density e when the matter filling the space is described by a viscous fluid energy-momentum 
tensor. We now match equation (|2ip rewritten as 

[ - | x (2 - x) + e - | Co x ] t- 2 + 

+ [±(y-x)(y-l) + ei - | Co^r e ~ 1/2 - §(j, - x) Co v^o~] t*" -1-2 = , (28) 

with the hydrodynamical ones T" v — 0. It is worth noting that, in the non-viscous case 
(Co = 0), the energy density solution is determined without exploiting the hydrodynamical 
equations, as in [2J, since e directly comes out from the 00-gravitational equation. In our ap- 
proximation (u a is neglected with respect to uq), the energy-momentum tensor conservation 
law provides the equation 

d t e + d t (\n^) [| e -Coe s 9 t (lnV7)] - 0, (29) 
which can be simplified as follows 

[2e (z-l)- fCozV^] t~ 3 + 

+ [e 1 (b{y-x-2)+2xb- lCoX 2 be~ 1/2 ) + 

+ §(y-a)&eo-fs(i/-aOCotyeo] *»---» = o . (30) 

Equations (|28p and (|30|) have to be combined together and solved order by order in the 
expansion in 1/t (in the asymptotic limit t — > 0). Since for the coherence of the solution, we 
impose y > x, by solving the leading-order identities we get 

X = i ' e « = | ^ ■ (31) 

1 4 C,0 

The parameter Co has here the restriction Co ^5 4 /3V3 in order to satisfy the condition x > 0. 
In this way the exponent of the metric power-law x runs from 1 (which corresponds to the 
non viscous limit Co = 0) to oo \1'S\ [2] • We remark that this constraint on Co arises from a 
zeroth-order analysis and defines the existence of a viscous Friedmann-like model, in which 
the early Universe has to expand with positive powers of time. 

Comparing now the two first-order identities (which involve the terms pro-portional to 
anc j we eas iiy get an algebraic equation for the y parameter 

y 2 -y(x + l) + 2x-2 = . (32) 

The solutions are y = 2, y = x — 10. Obviously the second one does not respect the condition 
y > x; hence the first order correction to the 3-metric is characterized by the following values 

y = 2, ei = -±x 3 + 2x 2 - 2x . (33) 

It is immediate to see that, in the non viscous case Co = 0, we obtain x = 1, eo = 
ei = — V 2 i which reproduce the energy density solution {SJ. 



1 We remark that in [2] (see Eq.(34), (35)) this solution is found by imposing the consistence of the a/3- 
Einstein equation and not as a pure dynamical condition derived by the solution of the perturbed hydro- 
dynamical equation. 
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By guaranteeing the consistence of the model, we now narrow the validity of the parameter 
x to the values which satisfy the constraint x < y. Thus, from (|3ip. the quasi-isotropic 
solution exists only if 

a <« - 55S . (") 

i.e., the viscosity is sufficiently small. For values of the viscous parameter Co that overcome 
the critical one (Co)> the quasi-isotropic expansion in the asymptotic limit as t — > can not 
be addressed, since perturbations would grow more rapidly than the zeroth-order terms. It is 
worth noting that the study of the perturbation dynamics in a pure isotropic picture yields 
a very similar asymptotic behavior when viscous effects are taken into account [EH [2] . The 
Fricdmann-singularity scheme is preserved only if we deal with limited values of the viscosity 
parameter, in particular we obtain the condition Qq S °^ < Co/3: this constraint is physically 
motivated if we consider, as it is, the Friedmann model as a particular case of the quasi- 
isotropic solution. 

Comments on the total pressure sign: The solution of the unperturbed dynamics gives 
rise to the expression of the metric exponent x in terms of the viscous parameter Co and to 
the zeroth-order expression of the energy density, which reads 

e =1J + -' ( 35 ) 

In order to characterize the effective expansion of the early Universe, let us now recall the 
expression of the total pressure p (|17p at leading-order: 

P = i e +TiCo\/ea;, (36) 

where we have used the 4-divergence (1201) truncated at zeroth-order. By using these identities, 
the condition p ^ yields the inequality 

Co < Co72 , (37) 

which strengths the constraint (|34|) and restricts the x-domain to [l, 4 /3]. 

The request to deal with a positive (at most zero) total pressure is consistent with the 
idea that bulk viscosity must not drastically change the standard dynamics of the isotropic 
Universe. In this respect, we address the domain Co ^ Co/2 as a physical restriction on the 
initial conditions for the existence of a well grounded quasi-isotropic solution. 

By concluding this section, we rewrite the expression of the energy density in order to 
analyze the density contrast evolution. In presence of bulk viscosity, e assumes the form 

3 a; 2 (x 3 /2 - 2x 2 + 2x) b 

r* ' (38) 

and, hence, the density contrast <5 can be written as 

S = -§ (x/4+ 1/x-l) bt 2 - x . (39) 

Since x runs from 1 to 2 as the viscosity increases towards its critical value, we note that, as 
in [inj , the density contrast evolution is strongly damped by the presence of dissipative effects 
which act on the perturbations. In this sense, we remark that bulk viscosity can damp the 
evolution of perturbations forward in time. This behavior implies that the density contrast 
approaches the singularity, i.e., 6 = 0, more weakly as t — > when the viscosity runs to Co- 
in correspondence with this threshold value the density contrast remains constant in time 
and hence it must be excluded by the possible Co choices. 
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7 The relation for the velocity and the 3-metric 



The 00-Einstcin equation provides the solution for the energy density; to perform a complete 
analysis of the quasi-isotropic model and to verify the consistence of our approximations, we 
now investigate the solutions of the Oa-components of the gravitational equations and the 
spatial a(3- ones. 

Imposing the condition s — 1 /2, the Einstein equation (I3bj) reads 

'- = t-«-w;»Jg + fe*V (40) 



Substituting (|38[) in the last equation, we get the following expression for the velocity, up 
to the leading-order of expansion (here in particular we neglect terms of order 0{t^ r ) and 

^ = ^(b a -^)t 3 -. (41) 

It is worth noting that, in our generalization, the assumption u 2 ~ 1 is well verified, since 
we immediately see that ii a u' J ~ t 6 ~ 3x , which can be neglected in the 4- velocity contraction 
u^u^ — —1; hence the approximated hydrodynamical equation (|30j) is still self-consistent 
using this expression of u a . 

Let us now write down equation (|3cp : here, the first two leading-orders of the right-hand 
side are 0(t~ 2 ) and (D(t~ x ) respectively only if x < 2 like in our scheme; hence u a u^ is 
neglected, as seen before, 0(t~ 2 ) terms cancel each other, while those proportional to t~ x 
give the following equation (which generalize (fT0| ) 



Pi + Abi + Bb5 f i + C<Sf = , (42) 
where the quantities A, B, C are defined as 

A=\(4-x 2 ), B=\{2x-l){x-2f-\x(x-2), C = -|(2 - x)(x - 1) , (43) 

respectively. Taking the trace of 1|42|). we obtain the relation (A + 3B) b = —P — 3C which 
provides the following equation 

2Ab^ p = (A + B)b, a , (44) 

when combined with the Ricci 3- tensor relation P^.g = 1 /2p-p. 

Therefore we are now able to write down the final form of the 3-velocity related to the 
perturbed metric tensor trace b : 

2 ~ X (A - B) t 3 - x b, a . (45) 



" \xA 

As it can be easily checked, the solution here constructed matches the non-viscous one (fTTj) 
if we set Co = and it is completely self-consistent up to the first two orders in time. As 
in the original analysis, the present model contains only three physically-arbitrary functions 
of the spatial coordinates, i.e., the six functions a a f3 minus three degrees of freedom ruled 
out by fixing suitable space coordinates. The only remaining free parameter of the model is 
viscous one, Co- 



8 Concluding remarks 

Our analysis outlined how the presence of bulk viscosity can deeply modify the quasi-isotropic 
solution in the asymptotic limit near cosmological singularity. The investigation started 
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from the modification of the Einstein equations, induced by a viscous matter term and then 
proceeded by the integration of the new gravitational equations matched together with the 
hydrodynamical ones, order by order in the 1/t expansion. 

As a main result, we have shown that the quasi- isotropic solution exists only for particular 
values of the bulk viscosity coefficient Co- When the dissipative effects become too relevant, 
we are not able to construct the solution following the line of the LK model. In fact, when 
Co approaches the threshold value Co = 2 /3\/3, the approximation scheme fails and the model 
becomes non self-consistent. 

By requiring that the viscosity parameter Co be under its critical value, we have also 
outlined how the behavior of the density contrast is deeply influenced by the presence of bulk 
viscosity. In fact, as far as dissipative effects are taken into account, the density contrast 
contraction (S — > as t — > 0), is damped until remaining constant if Co assumes its critical 
value. 

We conclude by stressing that our result is relevant near the singularity, where the volume 
of the quasi-isotropic Universe changes rapidly and as a consequence, the cosmological fluid 
has to follow this rapid variation by subsequent stages of thermal equilibrium. Then bulk 
viscosity emerges from the average non-equilibrium effects and it is expected to be increasingly 
relevant, when the singularity is approached. 
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